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AN INVESTIGATION OF AIRCRAFT HEATERS 
III - MEAST3RED AND PECEDICTED PERFORMANCE 
OF DOUBLE TUBE HEAT EXCHANGERS 
By R. C. Martiaelli, E. B. Weinberg, 
E, H. Morrill, and L. M, K, Boelter 

INTRODUCTION - 

Tvo double tube cylindrical heat exchangers, in which 
hot exhaust gases pass through the annular space and ven- 
tilating air passes through the center tube, have been 
tested to determine heat transfer performance and pressure 
drop. One of the exchangers was equipped with a smooth 
cylindrical air pipe (fig, 4), while the other utilized a 
dimpled type intensifier tube (fig. 5). 

The tests were performed in order 

1. To establish a simple, accurate method of pre- 

dicting the performance of the double tube 
heat exchanger, since this design is a basic 
element of practically all gas- air heaters. 

2. To compare the perfonuance of the straight and 

dimpled tubes. 

3. To determine the pressure drop across the units. 

EXPERIMENTAL EQUIPMENT 

In brief, the experimental equipment consisted of a 
source of hot gases and a source of cool air. These were 
metered and passed through the heater being tested. 

Measurement of appropriate temperatures and pressures 
allowed the calculation of the performance of the heater. 
Two sources of hot gases were utilized. The exhaust gases 
from a 50-horsepower gasoline engine were used at first, but 
it was discovered that hot gases from a natural gas burner 
yielded essentially the same results with much less opera- 
tional difficulty. The latter method is utilized at present. 
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Gasoline Engine Test Stand* 

An analysis of the thermal resistances existing in 
the double tube heat exchanger (see p, 7) indicated that 
the fundamental parameters governing the rate of heat 
transfer were the gas temperature and the weight rate of 
gas flow per unit area. The latter parameter may he 
varied by changing either the weight rate of gas flow or 
the cross-sectional area through which the gas moves. 

In absence of a large engine capable of producing 
high exhaust gas rates, a small (50 hp) engine was uti- 
lized to produce the hot gases, and high weight rates 
per unit area were obtained by the use of. a small exhaust 
gas annulus (0,D. = 3.07 in./ I.D. = 2.00' in,). 

The ventilating air was supplied to the heaters by a 
Eoots-type blower of 1000 cubic feet per minute capacity. 

The flow of exhaust gases was obtained by metering 
the fuel and combustion air to the engine. It was not 
necessary to meter the hot exhaust gases, since the 
weight rate per unit cross-sectipnal area is the important 
parameter which can be obtained by metering the fuel and , 
combustion air to the engine. The latter measurements 
were more easily obtained. The flow of ventilating air 
was measured by means of a calibrated lir-inch sharp-edge 
orifice. 

Temperatures were obtained at the following points: 

1. Ventilating air into heat exchanger 

2. Ventilating air out of heat exchanger (see p. 18) 

3. Exhaust gases into heat exchainger 

4. Exhaust gases out of heat exchanger . . 

5. Ventilating air tube wall (4 points) 

6. Outer tube of exhaust gas annulus (4 points) 

The exhaust gas temperatures were obtained by means- 
of the shielded thermocouples prqposed by the American 
Society of Mechanical Engineers (reference T). Chromel- 
alimel thermocouples were used throughout, 



Pressure drop along the ventllatixig air tube and the 
annular exhaust space were measvired "by means of U-tube 
vater manoineters,, Carefully installed piezometer rings 
were used on the air tube, but, owing to the construction 
of the gas annul us, single pressure tap holes were of ne- 
cessity utilised on the exhaust gas stream. 

Engine speed, load, exhaust gas analysis, barometric 
pressure^ and air humidity were also determined. The en- 
tire heater unit was surrounded by 2 inches of sand. 

A schematic diagram of the gasoline engine test stand 
is shown in figure 1, 



Natural Gas Test Stand 

The natural gas test stand (fig. 2) utilized a 
1,000,000 Btu per hour gas biirner to create the hot gases. 
The flame from the burner was conducted through a converg- 
ing section to the exhaust gas annul us. Part of the second- 
ary air was supplied by a small blower in order to insure 
complete combustion before the gases struck the venti- 
lating air tube. 

The gas flow was controlled by means of a centrifugal 
exhaust fan located^ downstream from the heat exchanger. 

The rate of exhaust gas flow was measured by means of 
a calibrated le- by 3-inch venturi meter placed between 
the heater and the exhaust fan. 

All other measurements were made in the same manner 
as with the gasoline engine test stand. 



SmSOLS 

Ag^ heat- transfer area of air side, ft 
Ag heat-transfer area of gas side, ft^ 

A,, surface area of outer tube of annulus, ft^ 
w . ' 

c-j unit heat capacity of air at constant pressure, 
Btu/lb ^ 

C constant 



t 
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D' diameter, ft. 

Dq^ inner diameter of air tube, ft . 

Dg hydraulic diameter of gas annul us, ft 

ei,e3 emlssivities of inner and outer surfaces of annulus, 
respectively 

f^ imlt theimal convectlve* conductance, Btu/hr ft^ ^ 

f^g^ unit thermal convectlve conductance on air side, 
Btu/hr ft^ 9F 

unit thermal convectlve conductance on gaa side, 
Btu/hr ft^ ^ 

f equivalent unit thermal conductance for radiation, 
Btu/hr ft^ 

shape modulus, the factor in the radiation equation 
vhich allows for the .geometrical position of the 
radiating surfaces 

Fj; emissivity modulus, the factor in the radiation equa- 
tion which allows for the non-Planckian character 
. of the radiation surfaces 

g gravitational force per unit mass, lh/(llD seo^/ft) 

G weight rate of flow per unit area, Ib/lir ft^ 

G^a weight rate of flow per unit area for air, Ib/hr ft^ 

Grg weight rate of flow per unit area for gas, Ib/hr ft^ 

kg^ thermal conductivity of air, Btu/hr ft^ (^/ft) 

N distance between pressure taps, ft 

P pressure, lb/ft ^ 

pressure at entrance to heat exchanger, lb/ft3 
pressur*e at exit from heat exchanger, lb/ft s 
predicted rate of heat transfer, Btu/hr 
q^ rate of heat transfer by convection, Btu/hr 
qj^ measured rate of heat transfer, Btu/hr 
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qp rate of heat transfer "by radiation, Btu/hr 

r , radial coordinate measured from center line of tube, ft 

r^ radius of pipe, ft 

tp average intensifier tube wall temperature, 

t^ average temperature of outer vail of annulus , ^ 

T absolute temperature, ^ 

Tg^^, mixed mean absolute temperature of air entering 
heating section, ^ 

Tq^3, Tg mixed mean absolute temperature of air leaving 
heating section, 

Ta arithmetic average absolute temperature of air in 
heater, 9R 

Tg arithmetic average absolute temperature of gas in 
heater, 9R 

Tp average absolute temperature of intensifier tube vail, 
T« average absolute temperature of outer wall of annulus, 

u velocity of fluid at any point r, ft/sec 

Uyn^ y velocity of fluid at center of pipe, ft/sec 

V mean velocity of fluid based on rate of discharge and 
pipe area, ft/ sec 

Vj^ mean velocity of fluid at entrance to heating section, 
ft/sec 

V3 mean velocity of fluid leaving heating section, ft/sec 

Wg^ air weight rate, Ib/hr 

Wg exhaust gas weight rate, Ib/hr 

X coordinate measiired along tube, ft 

y coordinate measured perpendicular to tube wall, ft 



y veight density of air at any temperature and pres- 
siire, lb/ft 2* 

7 veight density of air at. entrance to heating section, 
' lb/ft 3 

7 weight density of air leaving heating section, Ih/ft^ 

7rp weight density of air at temperature Ih/ft^ 

AP pressure drop across pipe, lb/ft ^ 

APm isothermal pressure drop due to friction at tempera- 
^ ture Ti, lb/ft2 

~ pressure drop per foot, (lb/ft^)/ft 
AL 

At difference in temperature between the tube wall and 
any point y in the fluid stream, ^ 



- (-^ga- -^aa)^ ^ 

At ^^y tube temperature minus temperature of fluid at cen- 
ter of stream, 

^"^mean temperature minus temperati;re of fluid after 

complete mixing, ^ 

At^j^ logarithmic mean temperature difference, ^ 

t isothermal friction factor 

^T^ isothermal friction factor for air at temperature T^^ 

^ Tg^ isothermal friction factor for air at temperature 

t rp isothermal friction factor for air at temperature T 

|i viscosity of air, lb sec/ft^ 

jjuj viscosity of air at temperature T, lb sec/fts 

Tg^^ mixed-mean temperature of air at entrance to heating 
section, op 

Tg^g mixed-mean temperature of air leaving heating section. 

Op 



mixed-mean temperature of gas at entrance to heating 
section, 9F 



IrrQ mixed -mean temperature of gas leaving heating sec- 
tion, ^ 

2 

Nu = . Nusselt modulus 

k 

Pj. ^ op g 5600^ Prandtl modulus 
k 

Re s — £^ , Reynolds modulus 

3600 |i g 

ANALYSIS OF MECHANISM OF HEAT TRANSFER 

IN DOXJELE TUBE HEAT EXCHANGERS 

In a double tube heat exchaiiger, in which the exhaust 
gas flovs in the annular space and the ventilating air in 
the central tube (fig, 3), the mechanism of heat transfer 
at any point along the length of the inner pipe proceeds 
as follows: 



tv 




Figure 3.- Cross section of double tube heat exchanger. 
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(a) Heat flows frcan the exhaust gas to the Inner tube 
by convection. The rate of convectlve transfer 
1b given by: 



die = fcg(Tg - t^)dA, 



(1) 



(t) Heat flows between the outer and the inner tube by 
radiation through a diathermanous* medium. 



dqtr = 0,173 



\ioo/ Vioo/' 



(2) 



For the tubular heat exchanger (reference 2, p. 54) 
■ F^ = 1.00 



F = 



— + -a { — 



(3) 



- 11 

/ 



An equivalent , unit conductance for radiation may be 
defined so that 



d(lr = fv,(T 



r^'g 



tp)dAg 



Thus, 



fr = 



0.173 Ip^ Fj; 



i2L 

.100 



y - m 



(4) 

(5) 



(c) Heat may l^e transferred to the inner tube by 

gaseous radiation from the H^O and COg pres- 
ent in the exhaust gas (reference 3, p. 297). 
A calculation reveals this quantity to be 
negligible in most cases, 

(d) All of the heat which is transferred to the outer 

surface of the inner tube (the sum of the quanti- 
ties discussed above) must be transferred to the 
air flowing through the inner tube by convection. 
Thus, if the gaseous radiation is neglected 



^Affording a free passage to radiant energy. 



(6) 



The three equations (l), (4), (6) allov the predic- 
tion of the thermal performance of the heater. For sim- 
plicity, however, it is customary to eliminate the tuhe 
wall temperature from the eqiiations (reference 4, p, XIV-1). 
Then, upon integration along the tube length* 

. ^ ±12 ^ , (7) 



r ^cg)^g 



^ca ^aj 



In eqLuation (7) 



(a) The logarithmic mean temperature difference may "be 
determined readily from the meas-ured magnitudes 
of the gas and air temperatures entering and 
leaving the exchanger. .Thus: 



At 



(T 



S2 



'as 



) 



(h) The areas Ag 
measurement . 



\^g2 - ^a2/ 



(8) 



In 

and A, 



are readily determined by 



(c) The radiant unit conductance may be calculated by 
means of equation ( 5 ) . 

The magnitudes of the unit convective conductances 
^ca> ^cg ^^3"*^ ^ow be predicted. 

For the turbulent flow inside a circular tube, by 
reasoning based on an idealized analogy between heat and 
momentum transfer, it may be shown (reference 5, p. 447) 
that : 

*The use of the log mean temperature difference AtjjQ, is 
based on the postulate that (fj. + fcg) and fca are 
independent of the tube length (as well as several other 
restrictions). (See reference 4, p..XrV-l.) Although all 
these postulates are not exactly satisfied by the experi- 
mental equipment, the utilization of the log mean temper- 
ature difference will be found to be sufficiently accurate. 
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Nu 



Ee Pr G c. 



5.0 



8 Atjngan. 
Pr.+2n(l + -5 Pr) + 0.50 In f~ 



1^0 V 8/ 



(9) 



where 



friction factor for smooth tuhes defined "by 
4P ^ ^ L v2 



D 2g 



At ^fly tuhe temperature - temperature of fluid at center 
of stream, ^ 

^"^mean ^^^^ temperature - temperature of fluid after com- 
plete mixing, ^ 

Pr Prandtl modulus 

Ee Reynolds modulus 

Nu Nusselt modulus 

A close empirical approximation of the more exact 
equation (9) was presented by Colburn (reference 6, p. 174) 



j£- Pr' . 1 (10) 

a Cp 8- 

Until further data are made available the friction 
factor to be utilized in equation (10) probably should be 
that for smooth pipe, even when analyzing fairly rough 
tubes, for the extra pressure drop produced by roughness 
does not appear to yield a proportional increase in fo. 
(See reference 7, p. 99,) 

The dimples in' the dimpled tube have an appreciable 
effect on the pressure drop, but only a small effect on 
the rate of heat transfer - for reasons other than the 
increased friction, as is discussed in a later paragraph. 
The question of the exact effect of roughness on the rate 
of heat transfer needs further clarification. 

For the purpose of analysis, the friction factor of 

Nikuradse (reference 8) was utilized. For a range of the 

Reynolds modulus of 10,000 to 100,000 the equation of 

Nikuradse may be closely approximated (reference 9) by 
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i =-2J^ (11) 

Solving for from equations (10 ) and (11) in terms 

of the physlcai properties of the gas**^ 

f^ := 0.948 X lO-4(c^O-333 |^--0.467 V t^V66?j G^^! (i2) 

Plotting the group {c^o ^:z:5 ^^o.4e7 k^o.ee?) ^ 

function of temperature and expressing the result as a 
power function of temperature allows the prediction of 
at high gas temperatures. This procedure yields 

fe = 5.56 X ICT^ T^'^^^^l^- (13) 

I ■• •• ^ • - 

Equation (13) may be utilized directly to determine 
^ca^ using the inside dimeter of the ventilating air 
tuhe for .. A. question arises, however,, concerning the 
value of D for the annular space. Recent investiga- 
tions (references 11 and 12) show that the hydraulic diam- 
eter,** that is, 4 X cross-sectio na l area of flow ^ ^j^^^^ 

. wetted perimeter 

utilized in the Reynolds number is a good criterion of the 
turbulence in the annulus. Substitution of this dimen- 
sion in equation (13) will be found to yield satisfactory 
results. 

Some recent data (references 12 and 15) indicate that 
an increase in fcg is to be expected, owing to the "annu- 
lus effect.** Because this correction is small 

*The properties of air (reference 10) were utilized for both 
the ventilating air and the exhaust gases. This procedure 
for exhaust gases was resorted to in absence of more precise 
data; but it is probably sufficiently accurate since the 
exhaust gas is ccaaposed mainly of nitrogen, and combustion 
is complete. 

**As contrasted, to the quantity 4x ^ross-sectional area 



heat-transfer perimeter 
which was suggested by Jordan (reference 13) and Nusselt 
(reference 14). 
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for the heater tested and there is a lack of agreement in 
the literature as to the exact value of the increase, this 
refinement was not included in -the' analysis. However, for 
exchangers in which the diameter of the outer siirface of 
the annulus is much greater than the diameter, of the air 
tube, this correction should not' be overlooked. 



Recapitulation . 

For prediction of the performance of simple gas-air, 
double tube heat exchangers the .following equations are 

recommended:. 



(a) (la .= -p 



(^r ^cg).^g ^ca -^a^ 
0.173 Fa Fe 



X 4 y X 4 



•^W ! -^p 

V 100 y ^ \ looy ^ 



(Tg- tp) 



(14) 



(15) 



(c) f 



,,ca 



(d) f, 



eg 



■ A 0.80 

5.56 X 10n^ T^o..9e^_ 



Q 0.80 

5..56-X. ler^ V r^o: 



(16) 



(17) 



The above equations were utilized to predict the rates 
of heat transfer for both the straight and dimpled tubes. 
This technique ^underestimates the (iimpled-tube .performaince 
slightly. This discrepancy is discussed in a later para- 
graph. •■ ' ■ " • 



SAMPLE CALCULATIONS 

Run no. F-V2 (Straight tube with mixing chamber) 
Rate of air flow, Wg^ = 382 Ib/hr 

Rate of gas flow, Wg-^= 379 Ib/hr 
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Inside diameter of air tube, Da - 0.149 ft 
Hydraulic diameter of gas annulus, Dg = 0.0892 ft 
Heat transfer area, air side, Aa = 2,23 ft^ 
Heat transfer area, gas side, Ag = 2,49 ft^ 
Cross- sectional area of air flov, 0,0175 ft^ 
Cross- sectional area, of gas flow, 0.0295 ft^ 
Air tomP'--:rtturc entering, ig^.^ = 130*^ F 
Air temperature (mixed mean) leaving, T^^g = 416<^ F 
Gas temperature, entering, T^,^ = 1626^ F 
Gas temperatiirej leaving, T^^ = 1243^ F 

(a) Calculation of the log mean temperature difference: 

At, = (1626 - 15C) - <:i24 5 - 416) ^ p 



(t) Calculation of f, 



ca 



Tg, = "^ai '^ag ^ ^ 7330 



W 

G„ = = 21,800 Ih/hr ft^ 

* 0.0175 • • 



Da = 0.149 ft 

p 0. 80 

f,a = 5.56XlO-X°*''' 
ca a 0.20 



= (5.56) X 10-^733)°-""^i21800)lff^ 

(0.149)°-^° 



=16.8 Btu/hr ft^ Of 
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(c) Calculation of f^g 

^ 2 
¥ 

G„ = = 12,800 iWlir 

0.0295 



Da = 0.0892 ft 



fcg = 5.56 X 10-^Tg°'=^^. 



=16.1 Btu/hr fta 



n 0.80 

8 



(d) Calculation of f^. 

The average temperature of the ventilating air tube 
(tp) is required in order to calculate fj.. 



and 



Thus 



la = *'ca Aa ("tp - 

la = ^a °p ("^aa - V ) 



■'ca ^a 



= 582 X 0.241 (416 . 130) + 278 
17.2 X 2.23 

= 962° F 

Experimental data indicated that the average temper- 
ature along the outer surface of the exhaust gas annulus 
(ty) was 250° F less than the arithmetic average gas tem- 
perature. Thus: 

t^ = (Tg - 250) = 1435 - 250 
= 1185° F 



The emissivity of oxidized steel is approximately 
(reference 2, p, 46) eg.ual to 0,79 



0.79 5.87 \0,79 J 



©1 



= 0.720 
Ip^ = 1.00 

Then 

0.173 X 0.720 



1185__+_460 
V "100 ' / 



962 + 460 
. 100 > 



1435 - 962 

8..30 

(e) Prediction of 



4a = 



(f J, + .fgg)Ag Ag^ 



1135: 



1. 



(8.30 + 16.1) 2.49 16.8 X 2.23 



1135 



0.01645 + 0.0267 



a 26,300 Btu/hr 
The oheerved rate of heat transfer Xbi-. 

« 582 K 0.241(416 - 130) 
= 26,300 Btu/hr . 
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The ratio for this run is 



q.« 26300 

— = ~ = 1.00 

26300 

DISCUSSION OF RESULTS 



The data and results are presented in tatle I. 

Inspection of table I reveals that the rates of heat 
transfer predicted by means of equations (14)^ (15)^ (16), 
and (17) check the measured magnitudes for the straight tube 
wi.th an average deviation of 6 percent. The average of the 
ratio of the predicted to the measwed rates of heat 
transfer is 0.99. It should be noted that the change from 
the gasoline engine test stand to the nat^iral gas. test stand 
difl not affect the accuracy of the results. 

The equations (14), (15),. (16), and (17) slightly 
underestimate the therml capacity of the dimpled tube. The 
average of the ratio of predicted to measured rates. of heat 
transfer for the dimpled tube is 0,925 vith a mean deviation 
of ±3 percent. 

It may be concluded from these results that the dimples 
Increase the capacity of the dimpled tube about 7 percent, 
compared with the straight tube. 

This increase is to be expected for several reasons. 
First, the dimples decrease the cross- sectional area of flow 
for the ventilating air. Equation (16) indicates an increase 
in fca owing to this change. Secondly, the dimples 
provide an Increased surface area for heat transfer. Lastly, 
the dimples may increase the radiant. heat transfer by acting 
as pai-tial Hohlraums (ideal absorber). ' 

Further inspection of table I shows that the radiant 
heat transfer from the outer wall of the gas annulus to the 
ventilating air tube was appreciable in the experiments. 
This condition, however, may not exist in actual installa- 
tions where the outer wall of the annulus is at a lower 
temperature than existed in the laboratory. Thus the rates 
of heat transfer obtained in the laboratory are somewhat 
larger than may be expected under flight conditions for the 
same rates per unit area. 
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In order to help visualize the data presented in taMe 
I, a curve of 



4, 



At 



2m 



a ,r ' 

' 1 ^ lr::. r 

(f^ + fcg)Ag fa ^^aj 



(U) 



Is shown plotted in figure 6 as a function of the ventilating 
air rate per unit area for the following^ average experimental 
conditions for thq straight tube: : ' 

Average gas temperature, Tg -= 1500<^ F = ^ r -' 

Average air temperature, T^^ = 230^ F 

Logarithmic temperature difference, ^"*^2m ^ ^^^0^ F 

Unit conductance for radiation, f^, = 7.0 Btu/hr ft^ 

Weight rate of gas flow per unit area 

Curve (1), Gg := 13,000 Ib/hr ft^ 

Curve (2), Gg :r 26,000 Ib/hr ft^ 

For comparison with the predicted curves, data for 
several eyperimental runs are shown plotted in figjurp 6, 

The data shown were corrected to a logarithmic mean 
temperature difference of 1050^ F by multiplying the 

measured q^^ by the ratio ( ^t^^^ order to allow 
direct comparison with the curves. 

The rates of heat transfer to be expected at an ex- 
haust gas rate of 26,000 Ib/hr ft^ are also shown. The 
increase is very small at low air rates, owing to. high 
thermal resistance on the air side of the central tube. 

A series of charts for the solution of equation (14) 
for rapid calculation of heater sizes (neglecting radia- 
tion) will be found at the end of the report. 
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TEMPEEATURE DISTBIBUTION IN AIR STREAM 



A very important point of experimental technique which 
has teen overlooked by nmerous investigators is the fact 
that, even at high magnitudes of the Reynolds modulus, a 
decided temperature gradient exists across an air stream 
"being heated in a tube. That this must be the case may be 
demonstrated by an analysis of the thermal resistances from 
the tube wall to the center of the tube. It may be shown 
that (reference 5) the temperature distribution in the 
turbulent core is: 

Pr + 2n(l + 5 Pr)+ 0,5 In |k /4" 
At _ oO r^ ^ b ^^gj 



Pr + Zn(l + 5 Pr) + 0.5 In ~ 



if 



Integration of the temperature distribution across the tube 
yields: 



^ 




Utilizing the equation for the velocity distribution pre- 
sented by Bakhmeteff. (reference. 16) gives the resulting 

ratio of ^^^P^^- as shown in figure 7, 
^*lnax 

Figure 7 reveals that, unless a very excellent mix- 
ing chamber is utilized, a single measurement of tempera- 
ture at the center of the air stream will introduce serious 
errors. Several runs were made purposely with no mixing 
chamber in order to demonstrate this point. 

The data are shown in table II. The apparent tem- 
perature of the heated ventilating air in these runs was 
based on the readings of a single couple at the center of 
the air stream with no mixing chamber. In contrast to the 
satisfactory results shown in table I, the omission of the 
mixing chamber yielded apparent rates of heat transfer 
which were on the average 30 percent below the predicted 
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values for the straight tube and 10 percent helow the pre- 
dicted values for the dimpled tube. The smaller discrepancy 
for the dimpled tube data is readily explained by a con- 
sideration of the temperature distribution vhich was measur- 
ed across the tube diameter for the straight and dimpled 
tubes (fig. 8) during the tests. An inspection of these 
curves reveals that the dimples act as a partial mixing 
chamber for the dimpled tube temperature distribution curve 
is much flat-Jier than that obtained in the straight tube. 
A comperiebn of the experimental curves and those predicted 
by equation (18) is shown. Application of equation (19) to 

the ciurves shown in figure 8 yields a ratio of ^.^^^^ = 6.80 

Atmax 

for the straight tube and 0.90 for the dimpled .tube, thus 
explaining the apparently greater temperature rise of the 
air passing through the dimpled tube. 

A direct comparison of the experimental results shown 
in table II would lead to the erroneous conclusion that the 
dimpled tube was 30 i>ercent more efficient than the straight 
tube, instead of the 7 percent improvement obtained by the 
more acctrrate mensuration of the exit air temperature. 



PRESSURE mOP 



The pressure drop along a tube in which a compressible 
fluid (gas) is being heated is not simply the frictional 
loss occurring along the tube, but also includes the effect 
of the expansion of the gas as it becomes heated. The lat- 
ter effect produces a large fraction of the total pressure 
drop, so that pressure losses observed under nonlsothermal 
conditions should be corrected to isothezmal conditions - 
or, preferably, isothermal pressiu:»e drops should be observed 
if the frictional loss through the heater is desired di- 
rectly. 

The analysis of the nonisothermal pressure drop through 
a cylindrical tube of constant diameter is presented below. 

The Bernoulli equation for a horizontal tube is: 
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For a constant weight rate there 'is a definite relation "bei- 
tween velocity and density. Thus, 



or 



Thus the Bernoulli equation l)epome8 : 
G dv " ^ 



J ■ 3600 g 3&00 2gD J 



By assumption of. a linear increase of velocity* vith length 
due to heating of th^ fluid, the velocity at entrance to 

/ T \ 

the heating, section is v^, -and at exit v^ ( =^ ) . Thus at 

.. A / 

any point x along the heating lerigth^ 

V = vi + vn^Ik - l\ i (23) 

where Tx and Ts are the absolute mixed mean tempera- 
tures of the fluid ;.entering and leaving the heating section. 

Substituting the above equation in the Bernoulli 
expression and integrating** yields:- 



Substituting 



3600 7i 

*The incr^ease in velocity is actually exponential, but as a 
first approximation a linear variation may be utilized. 

**Assuming t to be constant,. Actually ; varies slightly 
with lerigth due to the change in Reynolds modulus as the 
fluid becomes heated. 
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3600 72 



yields 

'*(3600)^2g ri D V^ax^ . 



(Pi - Pa) - 



(3600)^718 \Ti. 



(25) 



To obtain the iaothgrmal head ..loss at the temperature 
Ti, the friction factor t^' muat be replaced by the 

friction factor calculated a% til© temperature T, (Sm ) 
vThe; friction factor (reference 9) usually ±s expressed 



as: 



i 



T 



vO.2 0 



but (reference 10) 



Thus 



Tlius the filial equation becoimes: . , 



G8 



^ , 

'^1 (3600)^ 2g 71 V ■ 
(Pi - Pa) 



,1 .13 



(3600)^ 



/££ - l) 



(26) 
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where 

r T? ) = ^^T f isothermal pressure 

(3600)^ 2gr^ 1 drop due to friction 

at temperature T;^ 

A similar derivation is presented "by McAdama (reference 2-> 
p, 130). ... . • 

Thus the isothermal head loss due to friction^ based 

on the temperature Ti is equal to the measured non-. 
isothermal head loss (P^^ • less a correction 

l^'^s _ A and is multiplied by the ratio 

(3600)2 gVT^ J 

1,13' 

. Conversely, if the isothermal head loss is 




known, the nonisothermal loss may be calculated, • 

(Pi - Ps) = apt/ ^T""' + - ^) 



(27) 



It should be noted that the nonisothermal head loss is 
greater than the isothermal for a fluid heating, but less 
for a fluid cooling. 

Several runs were made in order to check the validity 
of the derivation shown above. 

The data and results are shown in tables III and IV, 
and are plotted in figures 9 and 10 i. A series of runs was 
made to obtain the isothermal pressure drops (table III), 
and another set of data obtained under nonisothermal 
(heating) conditions (table IV)* 

Inspection of figures 9 and 10 reveals that for both 
tubes the pressure drop**^ during the nonisothermal condi- 
tions is, rou^jhly, twice that during iBothormal conditions. 



*The pressure drop per foot is plotted for convenience, but 
it should be remembered that this value is only an average, 
for the pressure drop will not be a linear function of length 
due to heating of the ventilating air. 



23 



Application of equation (26) reveals that this discrepancy 
is due largely to the effect of the reduction of the air 
density upon heating, and the resulting increase in fluid 
velocity. Some of the liicreased pressure drop is due also 
to the increased air viscosity in the laminar sublayer 
(reference 5). 

Equation (26) vas utilized to predict the isothermal 
pressure drop from the nonisothernaal^ as follows: 

(a). The quantity, 

(Pi - I^) - ^ 




(ieod)' 



— 



was computed. This maQaitude is the isothermal pressure drop 
to he expected if the air passing through the heater is at 
the temperature T^, 

(h) From equation (26) it is noted that, if the Iso- 
thermal i^res sure drop at any other tOTiparature T Is de- 
sired,_ .'.*'"' 

^.(M^.f^'^''' (28) 

For direct comparison to the isothermal data which were 
ohtained at 78^ F all the magnitudes of APfj^ were correct- 
ed to 780 F hy means of equation (28). 

The isothermal pressure drops computed in this manner 
ftrom the nonisothermal. data checked the measured values with- 
in 27 percent for the straight tube and 7 percent for the 
dimpled tube. Thus, equation (26) dppears to be reasonably 
successful In correlating the isothermal and nonisothermal 
pressiure drops in tubular heat exchangers . 

As a further comparison, the Isothermal friction fac- 
tors for the dimpled and straight tubes are compared to the 
smooth pipe data of McAdsuns (reference 2, p. 172) and 
Nikuradae (reference 8) in flgiJtre 11. The magnitudes of the 
friction factor for the straight tube lie somewhat above the 
smooth pipe curve. The friction factor for the dimpled tube 
is about 1.6 times that for the smooth tube. This increase 
is to be .expected because of the effect of the dimples. 
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On the gas side of. the exchangers the foil oving aver- 
age nonlsothermal pressure drops were observed for a gas 
rate of 420 Ib/hr. 

Nonlsothermal Pressure Drop 

' ■ Ib/f t3 

ft '^81 82 

Straight tube 2.62 1500^ F . lOOOP F 

Dimpled tube 3v24 XSQOP F , 1000^ F 

The isothermal preseiore drop would be somewhat greater, 
as is shown by equation (26). 



CONCLUSIONS, 

1. The data presented in this section indicate that for 
testing the thennal output and the pressure drop performance 
of double tube heat exchangers, hot gases produced by a 
natural gas burner yield results substantially similar to 
engine exhaust gases . 

2. Eq^Tiations (14), (15), (16), and (17) are applicable 
to the predictions of the rates of heat transfer in double 
tube heat exchangers (exhaust gas to air), 

3. The dlmpled-type tube shows a slight advantage (7 
percent) in heat capacity over the straight tube. 

4. Unless a well-designed mixing chamber is utilized 
p.t the point where the temperature of a nonlsothermal gas 
streatn is measxired, appreciable errors will result In the 
measured temperature. ' 

•5. The nonlsothermal pressure drop along a tube in which 
a gas flows differs widely from the isothermal friction loss. 
Equation (26) allows the prediction of the isothermal loss 
from- the' nonlsothermal and equation (27) may be utilized to 
calculate the, nonlsothermal pressure drop when the Isotheimal 
is , known. 

• 6, 'The isothermal friction factor for the dimpled tube 
is 1.6 times that for the straight tube. The latter is some- 
what greater than the friction factor for smooth pipe. 

University of California, 
Berkeley, Calif. 
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APPEKDIX 

In order to facilitate the quick estimate of the heat- 
ir^ capacity of any double tube heat exchanger^ a series of 
cliarts have been dravn which allow the graphical determina- 
tion of heater capacities. These charts are based on equa- 
tions (14), (16), and (17). Eadiant heat transfer is 
omitted, . - 

Chart A allows the determination of the log mean 
temperature difference* At^^ X7hen the torminal ..temperature 

differences are known, = 

Chart B allows the deteminatidn of the thermal re- 
sistance of the heater per foot of lerigth 




when the. rates of gas flow, air flow, and 



the cross- sectional design of the heater are known. 

The thermal output of a given heater <Iq^ or the 

of heater required to give a; certain output N may be 
determined from the simple equation 




where 

N length of heat exchanger, ft 
P heat- transfer perimeter, ft. 

An example is presented to clarify the use of the charts. 
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Example 

Given a heater with the cross section shown in f ig. 12* 
The heater is 3 ft long. The following design conditions - 



are assumed: 

Rate of exhaust gas flow . . . ... • • • . 6000 Ib/hr 

Rate of ventilating air flow 2000 Ib/hr 

Temperature of exhaust gas entering 

heater 1600O F 

Temperature of exhaust gas leaving 

heater / 1100^ F . 

Temperature of air entering heater ^ -60O F 

Temperature of air leaving heater ...... 400^ F 

To determine: 



1. The output of the heater if the air flows in the 

spaces marked A in fig. 12, and the exhaust 
gases in the central space (B, C). 

2. The output of the heater if the flow is the same as 

under 1, but the central 3 -inch hole (C) is 
closed off. 

3. The output of the heater if the air flows in the 

spaces marked B and the exhaust gases through A. 

4. The average temperature of the heat transfer sur- 

face imder the above conditions. 

Notes: 1. Pressure drop and manifolding design considerations 
have been neglected. These, of course, may 
■ * invalidate some of the combinations proposed 

above. 

2. The assimed temperatiures will change with the 
heater output, but as a first approximation, 
the temperatures will be considered imchanged. 
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Solution - . '■ ' 

From f ig , ' XZ the f ollowing "data are olitalried : 



TMLE v.- HEiyrER DATA 





Proposal 

(1) 


Proposal 

(2) 


Proposal 
(4) 


CrosB-sectional area of air 
flow in ^ 


31.4 


31.4 


11,8 


Cross- sectional area of gas 
flow (sq. In. ) 


18.9 


11.8 


■31.4 


Wet t ed. per ime t er- of. ■ air -i 

flow (in.) 


105 


105 


89 


Wetted perimeter of gas 
flow (in.) 


80 


89 


105 


Heat-transfer perimeter for - 
air (in.) 


80 


80 


80 


Heat-transfer perimeter for 
gas (in.) 


80 


80 


80 



The average gas temperature is 1350^ F. The average air 
temperature is 170^ F/ The temperature., differences between 
gas and air are: 



At entrance Ati = 1660^ F 
At exit ■ 700O F 

From chart A: 

At iti = 1660^ F^ Atg = ^700^ Atj^^^ = 1100^ F 

From chart entering with the weight rate of air or 
exhaust gas flow, and proceeding as shown hy the dotted 
line^ the following values are obtained: 
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TABLE VI.- THERMAL RESISTANCES FROM CHAET B 









Proposal 

(3) ; 




• Proposal 
(1) 


Proposal 
(2) 




0.0175.. 


0.0175 . 


0,0062 •■, 






1 \ OF 


,0034 


.0020 


: . .0053 


V^Wgas (^Whr ftp 


C^c^/^alr Oc^ygas _ 




.0209 


.0195 


.0115 













Substituting in equation (Al) for a heater length Of 
N = 3 ft:* 




Proposal (1) qa = 157,000 Btu/lir 
Proposal (2) qa ^ 169,000 Btu/hr • 
Proposal (3) q^^- = 287^000 'Btu/hr 

Inspection of table VI reveals that when the' gas flows 
in the B spaces and the air In the A spades , the thermal 
resistance on the air side is so large that the heater ■ 
capacity is low. Closing off the central hole has little 
effect on the heater capacity. But, hy reversing the 
position of the air and the gas (proposal 3) the two thermal 
resistances are equalized and the heater capacity is practi- 
cally doubled. The optimum condition for heat transfer 
normally would he one in which the air and the gas side 
resistances are about equal. 



^Fov more accurate calculations the log mean temperature 
difference At^m should be recalculated to agree with the 
thermal output of the heater • 
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Chart B my be used to design for this condition. If 
the length of heater at capacity and the log mean tempera- 
ture difference are fixed, the total resistance 

may be determlnied. One-half of this 

total resistance will be the desired resistance on the gas 
or the air side. Utilizing chart by choosing various 
combinations of cross- sectional areas and perimeters, a 
heater may be evolved with the desired chaaracteristics. 

The average temperature, of the heat-transfer surface 
also may be estimated. By amlogy to electrical circuits 
(reference 17) the average temperature of the heat- transfer 
surface may be estimated from the resistances shown in 
table VI. 



Total resistance ... 0.0209 ^ 



Thus for proposal (1) 

Btu/hr ft 

Gas side resistance 0,0034 ~ 

Btu/hr ft 

Average gas temperature .... 1350^ F 

Average air temperature .... 170^ F 

Difference 1800^ F 

Thus the drop in temperatvire from the gas to the tube 
wall is 

8 P 0.0209 

or . tp := 1350 - 192 = 11580 F 

In a similar manner for 
Proposal (2) tp = 1229^ F 
Proposal (3) tp = 805^ F 
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Thus, the proper proportioning of the thermal resistances 
on the air and gas side not only improved the heater capacity, 
but lowered the temperatxires of the heat transfer surface 
over 400^ F. 

The foregoing" method, however, wiir not reveal local high 
temperature conditions which may exist in a heater tecause of 
improper design of entrance conditions, manifolding, etc. 

The atove example was selected to illustrate the use of 
the charts and Is not to he considered a suggested d-esign* 

The authors wish to express appreciation to Messrs. H. F. 
Brockscmidt, Trihus., D. Dubain, M, A. Miller, and F. Hamaker 
for their help in the performaiice of the tests and the prepa- 
ration of the report, and to Messrs. H. Eagles and H. Poeland 
for constructing the laboratory equipment. 
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TABLE III.- ISOTHERMAL PRESSURE DROP - AIR SIDE 



Run 


(llj/hr) 


( ) 
Vjxr ft*/ 




Measured 
pressure 
drop 


Heynolds 
numter 


Friction 
factor 

















Straight tute (distance 'between pressure taps, 6,16 ft) 



KIS-l 


3 94 


23.500 


78 


1.47 


75,400 


0.0269 


KIS-2 


411 


23 ,500 


78 


1.64 


78, 900 


.0272 


KIS-3 


561 


32 ,100 


78 


2.58 


107,500 


.0229 


KIS-4 


631 


36 ,100 


78 


3.14 


121,000 


.0221 


KIS-5 


681 


38,900 


78 


3.33 


130 ,000 


.0200 



Dimpled tute (distance between pressure taps, 6.21 ft) 



KID-l 


670 


34 ,600 


78 


4.42 


122,000 


0.0354 


KID-2 


617 


32 ,200 


78 


3. 64 


113 ,000 


.0354 


KID-3 


525 


27 ,200 


78 


2. 93 


96 ,000 


.03 82 


KID-.4 


456 


23 ,600 


97 


2.48 


81,100 


.0412 


KID~5 


408 


21 ,200 


97 


2.18 


73,000 


,0448 
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